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Abstract 

We have used the latest data from the ALEPH Collaboration to extract values 
for QCD condensates up to dimension d = 12 in the V — A channel and up to 
dimension d = 8 in the V, A and V + A channels. Performing 2- and 3-parameter 
fits, we obtain new results for the correlations of condensates. The results are 
consistent among themselves and agree with most of the previous results found 
in the literature. 
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1 Introduction 



QCD is widely considered to be a good candidate for a theory of the strong interactions. 
Asymptotic freedom allows us to perform a perturbative treatment of strong interactions 
at short distances. Long distance behaviour, however, is not fully understood: it is 
commonly believed that, due to the nontrivial structure of the physical vacuum, the 
perturbation expansion does not completely define the theory and that one has to 
add non-perturbative effects as well. In order to make a comparison with experiments 
possible even in the resonance energy range, Shifman, Vainshtein and Zakharov pQ have 
proposed to use the Operator Product Expansion (OPE) and to introduce the vacuum 
expectation values of operators occurring in the OPE, the so called condensates, as 
phenomenological parameters. It is worth to study these parameters in order to see 
whether one can indeed obtain a consistent description of the low energy hadronic 
physics and get more insight into the properties of the QCD vacuum. It is, in particular, 
important to determine the range of values of the condensates allowed by presently 
available experimental data. 

Condensates are needed in the description of two-point functions of hadronic cur- 
rents: together with the results of perturbative QCD they allow us to obtain approx- 
imate theoretical predictions for the hadronic current correlators in the space-like re- 
gion. On the other hand, in the time-like region, the discontinuity of these amplitudes 
is directly related to measurable quantities. Analyticity strongly correlates the energy 
dependence of the amplitudes in these two regions; however, due to errors affecting 
both the theoretical predictions as well as the experimental data, the relation of the 
two-point functions in the time-like and space-like regions must be carefully analysed 
before it can be used in applications. 

There are several methods, generically called QCD sum rules [HE1E], that have been 
used in the past for obtaining values of the QCD condensates. They all rely implicitly 
on the assumption that an analytic extrapolation from the time-like to the space-like 
region is possible without introducing additional uncertainties. Therefore they can 
include theoretical errors in the space-like region only at a qualitative level, and/or need 
(explicit and implicit) assumptions on the derivatives of the amplitudes. A quantitative 
estimate of the errors including both experimental and theoretical ones (truncation of 
the perturbative and operator product expansions) is therefore very difficult in these 
approaches. The application of fully controlled analytic extrapolation techniques should 
remedy these effects. There are a few methods of this sort, in which the error channels 
in the space-like region are defined through L 2 -norms [H E] or L°°-norms [El El E] • 

The functional method that we have developed and used in a previous publication for 
an analysis of the correlator of the V — A current [I] , allows us to extract within rather 
general assumptions the condensates from a comparison of the time-like experimental 
data with the asymptotic space-like results from theory. We will see that the price to 
be paid for the generality of assumptions is relatively large errors in the values of the 
extracted parameters. In this respect, our method is not superior to other approaches; 
however, we hope that our results provide additional confidence in the numerical results 
obtained with the help of methods based on QCD sum rules [H [21 El El fTUl [TTj . 
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The paper is organised as follows: In the first section the theoretical aspects of 
hadronic r-decays are considered while in section 3 we describe the experimental data 
at hand. The actual method used to extract the condensates is presented in section 4. 
Basically the method is the one reported in [I] but here a generalisation to all channels 
is performed. The results are summarised in section 5. We quote results in all the 
four channels where data are available. Results for the V — A channel were already 
published in [I] ; for completeness, we repeat those previous results in the present paper 
and extend our analysis by a new 3-parameter fit. New results were obtained in the V, 
A and V + A channels from 1-, 2- and 3-parameter fits. We shall also discuss important 
consistency checks of our method in section 6. The comparison of our findings with 
others present in the literature is performed in section 7. 



2 Theoretical description of hadronic r-decays 

The r lepton is heavy enough (m T = 1.777 GeV) to decay not only into other leptons, 
but into final states involving hadrons as well. These decays offer an ideal laboratory 
for the study of strong interactions, including the transition from the perturbative to 
the non-perturbative regime of QCD in the simplest possible reaction. This might 
explain the tremendous efforts ongoing in both theoretical and experimental studies of 
r physics. 

We consider the correlator of hadronic vector and axial-vector charged currents, 
Jn = Vp = vrtyd and J^ = A^ = uj^ 5 d, 

K = if d A xe^(TJ,{x)JM) (1) 

= (-g^ 2 + %<i») ni 0+1) (g 2 ) + M 2 4 0) (q 2 ) ■ 

The conservation of the vector current implies Ily = 0. The connection to experimental 
observables is most easily expressed with the help of the spectral functions which are 
related to the absorptive part of the correlators. Using the normalisation as defined in 
most of the previous publications, the functions 

Vj(s) = 27rlmn^(s), a^s) = 27rImIT^ ) (s) (2) 

can be extracted from the decay spectrum of hadronic r-decays. 

The hadronic polarisation tensor can be rewritten using the OPE: 

n (j),v,A 

d>0 V ' 

where Od = Cd(Od) is the short hand notation for the QCD non-perturbative condensate 
(Od) of dimension d and its associated perturbative Wilson coefficient Cd\ s = q 2 is the 
squared momentum transfer. 
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The contribution to ([3]) of lowest dimension, d — 0, is entirely given by perturbation 
theory. For that reason it is useful to separate the two contributions in ([3]): 



^A(s)=^iyJs)+Ug Ey J S ). (4) 

As will become clear when we describe our method (sectionHJ) it is enough to consider 
the first derivative of the perturbative part, i.e., the Adler function 

D(s) = - s ^n PT ( S ), (5) 

which is known in the massless-quark limit up to terms of order af. After re-summing 
the leading logarithms it reads, for space-like momenta (s < 0): 

^vM-l±Kj^X. (6) 

n>Q V 7 



The coefficients K n are the same for both V and A channels. For 3 flavours, in MS 
regularisation, K = K x = 1, K 2 = 1.64 [I2j EH E] , K 3 = 6.37 $M EE] and for K 4 
there are two estimates K± = 25 ± 25 [17J and K 4 = 27 ± 16 [18J. 

For the correlators with spin + 1 and in the chiral limit, the non-perturbative part 
has the form 

n8&» = E (^s (1 + cr°' v ' A ^) . (7) 

where perturbative corrections of order 0(a s ) are taken into account, described by 
coefficients c^ L °. Some of the NLO coefficients were calculated in [T9 l [201 f2Tj . 

The parameters Od can be related to vacuum expectation values of products of 
quark and gluon field operators [Ij. Often vacuum dominance or the factorization 
approximation, which holds, e.g., in the large- iV c limit, is assumed. Our analysis does 
not rely on such a representation. 



3 Experimental data 

Since its discovery, the r lepton has been studied with ever-increasing precision at ev- 
ery new e + e~ collider that has gone into operation. We are particularly interested in 
the comprehensive measurements of exclusive hadronic branching ratios from ALEPH 
[221 123] and of the non-strange spectral functions from ALEPH [22J [23] and OPAL 
[24J that have yielded important contributions to studies of perturbative QCD at low 
energies and, in particular, to the measurement of a s (m^.). Recent measurements of a 
set of semi-exclusive branching ratios by DELPHI are also available [25], but have not 
been analyzed to a similar extent. A number of exclusive branching ratio measurements 
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from BaBar have been reported already [26J and with more work invested in the un- 
derstanding of these higher-multiplicity final states we may expect high-precision data 
for the spectral functions to come also from the B-factories. 

We have chosen to use the final data from the ALEPH collaboration [23] because, 
as compared to those available from OPAL [21], they have the smallest experimental 
errors and provide a larger number of bins. 

The spectral functions ([2]) are obtained by dividing the normalised invariant mass- 
squared distribution of hadronic r decays, dR T y^/ 'ds, for a given hadronic mass \fs by 
the appropriate kinematic factor. They are then normalised to the branching fraction 
of the massless leptonic, i.e. electron, channel B e = (17.810 ± 0.039)% 



mz 



ViiS 



dR T . 



v 



Q\V ud \ 2 S Fj w B P ds 



mi 



mi 



ai(s) 



mz 



dR T 
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mz 
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mz 
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ml 



dR T 



6|Kd| 2 »S'EW B e ds 



1 — - 

mi 



-2 



Here Sew — 1.0198±0.0006 accounts for short distance electroweak radiative corrections 
[27] and the CKM mixing matrix element has the value \V ud \ = 0.9746 ± 0.0006 [21]. 
Due to the conservation of the vector current, there is no j = contribution to the 
vector spectral function, while the only contribution to a is assumed to come from 
the pion pole. The spectral function is connected, via partial conservation of the 
axial- vector current (PCAC), to the pion decay constant = 0.1307 GeV [29] through 
a^{s) = 2^pj{s-ml). 



4 The determination of condensates: a functional 
method 

Let us consider a set of functions F(s) which are admissible as a representation of the 
true correlator if they are real analytic functions in the complex s-plane cut along the 
time-like interval [so, oo) with sq > 0. The asymptotic behaviour of F(s) is restricted 
by fixing the number of subtractions in the dispersion relation between F(s) and its 
imaginary part f(s) = IhlF(s + iO) along the cut: 

, 1 r° f(z)dz 

F(s) = - / i±J — + subtractions . (9) 

Ws Z-S 

For our purpose it is convenient to get rid of the subtraction terms by taking an 
appropriate number of derivatives with respect s. We denote by C n (s,z) the kernel 
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occurring in the dispersion relation for the n-th derivative of F(s): 
1 f 00 

-s n F^\s) = - C n (s,z)f(z)dz. (10) 

For example, in the special case of the V — A correlator which vanishes identically in the 
chiral limit to all orders in QCD perturbation theory there are no subtractions needed 
and thus one takes n = 0. In this case, the dispersion relation will be identical to the one 
of Eq. © with no subtraction terms. On the other hand, the V, A and V+A correlators 
are dominated by their perturbative contributions and there is one subtraction needed 
in the dispersion relation. To get rid of this usually unknown constant term one needs 
to take the first derivative and set n — 1. 

In order to determine F(s) and f(s) we use the following two available sources of 
information: 



experimental data measured in the time-like interval r exp = [s , s max ]: 

, / \ 1 / v i( s )i m the V channel, , . 

fcxp[S) ~2^\ ai{s) + a {s), in the A channel. [ } 

The extension to V + A channels is straightforward, 
theoretical model given by perturbative QCD, i.e., 

— the prediction for F(s) in the space-like interval Tl = [s 2 , si]: 

F Q cd(s) = n^ 1 ^), J = V, A, V + A, V -A (12) 



- and /qcd(s) = ImF QCD (s + «0)| sg(Smax)Oo) since QCD is expected to be reli- 
able for large energies. 

As a next step in extracting values for the condensates, we split the integral on the 
r.h.s. of the modified dispersion relation ffTUj) into two parts: one that can be described 
by the experiment and the other one by the theoretical model, i.e., QCD: 

-^qcd(s) - - / Cn(s, z)f qC n(z)dz =- C n {s, z)f(z)dz . (13) 



QCD prediction: Fq CD (s) experimental data 



The goal of the method is to check if there exists a function F(s) which is in accord 
with both the data on r exp and the model on i\. For doing this, one can use an L 2 - 
norm approach and define two functionals xlLf] an d XrII]- XrU] compares the true 
amplitude f(s) with the data. Here one can take into account not only experimental 
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errors on each individual bin, but the full correlation of available data by using the 
covariance matrix V provided by ALEPH as a weight function. Therefore we define 



xlif] 



exp | 



dx 



dx'V-\x, X')(f(x) - f exp (x))(f( X ') - /exp(x')). (14) 



As a measure for the agreement of the true function f(s) with the theory, we define 
Xz,Lf] by comparing the left and right hand sides of (TT3T) 



xlif) 



w L (x) 



-^qcd( x ) 



7T 



C n (x, x')f{x')dx' I c/x 



so 



(15) 



where is a weight function for the space-like interval, i.e. an a-priori estimate of the 
accuracy of the QCD predictions, and written as l/a 2 L (s). should be chosen as a 

continuous, strictly positive function of s G Tl and encodes errors due to the truncation 
of the perturbative series and the OPE. It is expected to decrease as \s\ — > oo and diverge 
for s — > 0. For example, in the case of the V — A correlator we will use the next higher 
dimension contribution in the OPE as an error estimate on the space-like region. In the 
case of the V, A and V + A correlators the situation is a bit more complicated since they 
are dominated by their perturbative part. Thus one has three possibilities to define an 
error corridor: one can use the last known term of the perturbation series, or the first 
omitted term in the OPE, or a combination of the two of them. As an illustration, in 



the case of a 1-parameter fit a L 
1 



V,A, 
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cr, 
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(16) 
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(-x) 3 



The factor 3 in front of 0^' A in the expression of a£ A (x) comes from the fact that in 
the V, A and V + A channels one needs to set n — 1 in Eq. ffTUl) and thus take the first 
derivative of the operator product expansion in Eq. (GO). K 3 was given in section El 

In order to find the true function /(s) one can combine the information contained 
in these two functionals by means of Lagrange multipliers and find the unrestricted 
minimum of 



subject to the condition 

Xlif] < Xexp = ^V(«i,*)V'(si,a;)V- 1 («i, s 



(17) 



1,3 
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which will be the criterion to determine the Lagrange multiplier /i. This procedure 
leads to an integral equation for the imaginary part of the true amplitude, f(x; fi): 



f{x\n) = / e xp(z) + 




(19) 



+A 




dzK(x,z)f(z;ii), 



where A = l//i and 



JC(x, z) 



7T 





dyV(y,x) / dx' w L (x')C n (x' , y)C n (x', z) 



(20) 



Thus, the algorithm for determining acceptable values for the condensates is the 
following [30] : 

• Choose a model by stating how many terms in the OPE should be taken into 
account. The term with the next-highest dimension is used to define an error 
corridor in the space-like region; 

• Solve the integral equation (JT9l) iteratively, with fixed values of the free parameters 
of the chosen model, until the Lagrange multiplier fi satisfies the condition ffT8l) ; 

• Calculate xl corresponding to the above solution as a function of the free param- 



• Minimise Xl with respect to variations of the model parameters. The correspond- 
ing parameter values are the condensates we are looking for; 

• Determine confidence regions around the fitted parameters by solving 



Here we will assume that the underlying probability distribution is Gaussian, 
fixing Ax 2 to reflect the conventional 1-, 2- and 3<r-contours for n-parameter fits. 
Numerical values for errors will be given for la confidence intervals. 

5 Results and discussion 

The most important step of our analysis is to find a function /(s) (Eq. (|T9"]) ) which 
provides a best fit to both the data and the theoretical model. A direct comparison of 
the experimental data with the regularised function f(s) obtained from 1-parameter fits 
is shown in Fig. [TJ We find nice agreement over the full range of s with the exception of 
the highest s-bins. Here the spread of data points is apparently wider than individual 
errors on single data points. The largest differences are visible in the V + A channel 



eters; 



(21) 



S 



at s > 2 GeV 2 where the discrepancies in the V and A channels between data and 
the fitted function f(s) get enhanced while they appear to be compensated in the 
V — A channel. We emphasize that we have used the full correlation matrix provided 
by ALEPH. Correlations are certainly important in our fit, however, they can not be 
visualized in our figure. 
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Figure 1: The regularised function f(s) (see Eq. ( IT9"1) ) compared to the 
experimental data of Ref. [22]. f(s) was obtained with 1-parameter fits using 
the 120 first s-bins, i.e. from s < 3 GeV 2 . 
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5.1 1-parameter fits 

Let us start with 1-parameter fits and quote results for condensates of dimension d — 4 
(V, A, V + A) and d = 6 (V - A) @]: 

0\ = 1.6±g;t x l(T 3 GeV 4 , x!, min = 49, 

Of = 2.6+^ x l(T 3 GeV 4 , xi, min = 2.8, 

0™ = 4.218| x KT 3 GeV 4 , xi, min = 19, 

0\~ A = -5.9±i;S x 10- 3 GeV 6 , ^ min = 0.17. 

In the 1-parameter fits we have fixed all higher- dimension condensates to be zero. The 
results of the V — A analysis had been given earlier [4]. There we obtained an acceptable 
fit when choosing an error corridor defined by the dimension d = 8 condensate in the 
space-like region with 

\OY~ A \ ~ 1.3 x 10~ 3 GeV 8 . (22) 

I ° I max v ' 

Our fit thus provides an indirect estimate of the upper limit of \Oj |. In the other 
channels we have used the coefficient K 3 in the perturbative expansion of the Adler 
function to define the theory error. 

5.2 2-parameter fits 

When performing 2-parameter fits we aim to determine simultaneously the first two 
relevant condensates appearing in the operator product expansion. The additional 
freedom in the fit provided by the second parameter allows us, in general, to obtain 
better fits. In fact, it turns out that, except for the V — A channel, the condensates of 
next-to-lowest dimension have significant non-zero values, in contrast to the assumption 
underlying the 1-parameter fits. 

In Fig. [2] we show contours of constant xt- O ne can see that we find strong correla- 
tions between the two free parameters. This correlation allows us to determine a linear 
combination of the two parameters with a well defined and rather small error: 

0\ + 0.65 GeV 2 0\ = 0.66j£! x 10- 3 GeV 6 , 
O a + 0.65 GeV 2 O a = 1.60l°l x !0~ 3 GeV 6 , 

(23) 

OX +A + 0.65 GeV 2 O v A +A = 2.20l°f x lO" 3 GeV 6 , 

0%~ A + 2.22 GeV 2 0^ A = -18.30l°;i| x lO" 3 GeV 8 . 

The location of the minima, i.e. the central values of the fitted parameters are quoted 
in Tab. [U As expected, the value of 0^~ A found in the 2-parameter fit has the same 
order of magnitude as the corresponding estimate found from the 1-parameter fit (see 
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Figure 2: 2-parameter fits: 1-, 2- and 3a confidence regions defined by 
contours of constant x\ ' n tne P' ane of the two fitted parameters. The central 
values (see Tab. 15. 2D are marked by dashed lines. 



Eq. 1221) . Similarly to the 1-parameter fit, we can now give an estimated upper limit of 
the V — A condensate of dimension d = 10 which was used in the 2-parameter fit to 
define the error channel: 

\OY - A \ ~ 5.7 x l(T 3 GeV 10 . (24) 

I iu I max v ' 

As before, the error corridors for the V, A and V + A channels were defined by the 
perturbative contribution of order (Eq. [TBI) . 

Despite of the poor agreement of theory with data in the V and V + A channels, 
reflected by the large Xl mm values, it is important to remark that a consistent over-all 
picture has emerged from our fits. The results of the 2-parameter fits are in agreement 
with those from the 1-parameter fit. However, since here the values of 0^' A,V+A and 
0\~ A were left unconstrained, we have found larger ranges for 0^' A ' V+A and 0\~ A . 
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d = 6 


o o+O.T 
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-6.81*;° 


d = 8 








o 9+2.8 

O.Z_g 2 


X,L,ram 


20.4 


0.47 


7.1 


0.37 



Table 1: 2-parameter fits: central values of the fitted parameters in units of 
10~ 3 GeV d and the corresponding values of x| min . In the V, A and V + A 
channels the fitted parameters were the condensates of dimension d = 4 and 
6 while in the V — A channel we have fitted the dimension d = 6 and 8 
condensates. 

Note in particular that the slope of the correlations in the V, A and V + A channels is 
the same. For the V + A case, the values for both condensates of dimension d = 4 and 
d = 6 agree with the values found by actually taking the sum of the results from the V 
and A channels. 



5.3 3- parameter fits 

A 3-parameter fit is also possible, where we consider as free parameters the first three 
relevant condensates in the OPE. We have chosen to display our results as x|-contours 
in the planes defined by the three possible pairs of fit parameters. These correlations 
are shown in Fig. [3] for the V and A channels and in Fig. H] for the V + A channels. 
In every case we display 2-dimensional slices of the 3-dimensional allowed parameter 
ranges keeping the third parameter at its central value as obtained from the 3-parameter 
fit. These central values are summarised in Tab. El 

An estimated upper limit of the dimension d = 12 condensate in the V — A channel, 
the one used to define the error corridor, is 

\OYf A \ ~ 47 x 10" 3 GeV 12 (25) 

I lz I max v ' 

which is the expected order of magnitude. 

Again we observe that 3-parameter fits turn out to provide better x\ min va l ues than 
the 2-parameter fits. Obviously, the improved results are obtained since the higher- 
dimension condensate can be chosen non-zero in the fit. As a consequence, the central 
values of all condensates are shifted. In addition we observe that the 3-dimensional 
contours are not always ellipsoids and non-Gaussian errors play a role. 
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Figure 3: 3-parameter fits in the V and A channels: 1-, 2- and 3a confidence 
regions defined by contours of constant xi ' n tne plane of two of the fitted 
parameters and located at the central value of the third one. The central 
values (see Tab. [2]) are marked by dashed lines. 
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Figure 4: 3-parameter fits in the V — A and V + A channels: 1-, 2- and 3a 
confidence regions defined by contours of constant xi ' n tne plane of two of 
the fitted parameters and located at the central value of the third one. The 
central values (see Tab.|5]) are marked by dashed lines. 
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Table 2: 3-parameter fits: central values of the fitted parameters in units of 
10~ 3 GeV d and their corresponding values of x| min . In the V, A and V + A 
channels the fitted parameters were the condensates of dimension d = 4, 6 
and 8, while in the V — A channel we have fitted the dimension d = 6, 8 and 10 
condensates. The error estimates are obtained by projecting the 3-dimensional 
ranges allowed by the fit onto the corresponding selected parameter. 



6 Consistency checks 

In the following we present additional details of our algorithm and describe a number 
of consistency checks. We have, in particular, studied the behaviour of the algorithm 
and its results with respect to variations of parameters appearing in the analysis: the 
number of experimental data points N used for the extraction of condensates, the end- 
points of the time-like interval Tl, s\ and S2, as well as the dependence on the error 
parameter needed to define <jl- In this section, we restrict ourselves to the case of 
1-parameter fits. 



10 



O 



x 6 
> 

<D r, 
^ 2 

o 




-10 



Figure 5: Dependence on the number of data points. N denotes the number 
of the highest data bin used in the analysis. 
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Figure 6: Dependence on the lower end, s 2 , of the space-like interval: the 
fitted parameter as a function of s 2 (left) and \\ mm as 3 function of s 2 (right). 



One can show that the information on the condensates is contained in the lower part 
of the spectrum by adding or removing data points at largest s. In Fig. [5] (left panel) 
we show how the fit result for 0\~ A depends on the number of data points. One can 
observe a fast stabilisation of the result already for the N ~ 100 lowest-s data points. 
In contrast, for the case of the A-correlator, cf. Fig. [5] (right panel), one can observe 
that including or excluding data points above N = 120 has a stronger effect on the 
condensate Of. In this region the rapid oscillation of the data points as well as large 
experimental errors play an important role. The decision not to include experimental 
results from the highest bins in the analysis, is supported by inspecting Fig. [TJ there 
we found that the regularised function obtained in our analysis does not describe the 
data in the large-s region. We have thus decided to cut off data points above s = 3 
GeV 2 , i.e., we use only the N = 120 first data points. 

Figure [6] shows the behaviour of the algorithm with respect to changes of s 2 , the 
lower limit of the space- like interval T L . One should choose \s 2 \ as large as possible, but 
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we observe stability of the algorithm for values certainly not larger than a few times 10 2 
GeV 2 for the V — A analysis and even smaller values are required for the fits of the A 
condensate. This is illustrated in the right column of Fig. [U]which shows the dependence 
of xi mm 011 s 2- One can observe a plateau for the values of 0\~ A and Of as a function 
of S2 and thus infer the values used in the analysis to be s 2 = — 1 50 GeV 2 for V — A and 
s 2 = —3.5 GeV 2 for the A channel. For larger values of | 1 the minimum of \\ becomes 
larger than 1, signaling a bad fit. This behaviour may be due to numerical instabilities 
and limitations of experimental data; more important, however, is the fact that in 
our LO analysis we did not take into account perturbative higher-order corrections: 
the large perturbative tails of the V- and A-correlators become increasingly important 
when increasing \s2\ and the sensitivity to the low-energy condensates is lost. 

In the case of the A-correlator, when studying the dependence on s 2; we found that 
the best simultaneous description of experimental data and theory is obtained when 
we choose to define the error corridor with the help of the last known term in the 
perturbation series, i.e., 



Or \X) 



4tt 2 



O.s 



-X) 



71 



(26) 



In contrast, for an error corridor calculated from the next-higher term in the OPE, 
i.e., using a A (s) = const/s d , we observe a less distinct plateau when varying $2 and no 
stability for the results for Of. The fact that with the choice f[2"6"]) we find very stable fit 
results for Of even when increasing \s 2 \ beyond the range of the x\ m i„-plateau, makes 
us confident that our results for the A-correlator are meaningful. 

It is, in fact, to be expected that a definition of the error corridor with the help of 
the higher- dimensional terms in the OPE would narrow too fast (with a power of 1/s) 
and not leave space enough for perturbative contributions that fall only logarithmically. 
In contrast, for the V — A channel where perturbative contributions are absent, it was 




Figure 7: Dependence on si, the upper limit of the space-like interval T L for 

ol~ A . 
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possible to choose tighter error channels given by the omitted next-higher OPE term. 
We note that the definition of the error channels was the same for all 1-, 2-, and 3- 
parameter fits in the case of the V, A and V+A channels, whereas for the V — A channel 
we had to re-define the error channel when increasing the number of free parameters. 
This explains the observation that xi mm * s n °t necessarily increasing when going from 
1- to 2- and to 3-parameter fits for V — A condensates. 

There exists also a well defined plateau for the fitted parameters as a function of 
si. In the analysis we have chosen to use the values si = —1.0 GeV 2 for V — A and 
si = -0.4 GeV 2 for the A channel (see Fig. [7J). 

All these consistency checks were performed for the V and V + A channels as well 
with similar results and we found no justification to change the number of data points 
N used in the analysis or s 1; the upper limit of the space-like interval Fz,. Also, the 
dependence on the lower limit S2 has shown that the best simultaneous description 
of theory and data corresponds to an error corridor defined by the last known term 
in the perturbative series. For the correlator of the vector current, however, we are 
not convinced that we have obtained trustworthy results: first, ximin * s l ar g e even 
for a 3-parameter fit and, second, the fit results for 0\ and 0\ change by more than 
the estimated uncertainties when changing the number of free parameters in the fit. 
Therefore we do not discuss results for the V^-correlator further. Note, however, that 
0\ and Of are predicted to be equal, and fit results for one can be used to determine 
the other. 



7 Comparison with other results and conclusions 

There exists a number of previous extractions of QCD condensates in the literature, 
mainly based on sum rule approaches. For the V — A channel they are listed in Tab. [3] 
together with a repeated collection of our results. 

In most cases, errors given by the authors are in the order of 25%, sometimes even 
as small as 10%. For the d = 6 condensate, our results fall nicely in the same range, 
also with an error estimate which is comparable to that of other analyses. The spread 
of the central values is, however, larger than the typical error. We believe that the 
observed variation of these results represent the ambiguities inherent in the QCD sum 
rule approach. 

The situation is more difficult to summarize in the case of the higher- dimensional 
V — A-condensates: the variation of results from different analyses is even bigger, 
but estimates of relative errors are again in some cases similar to those of the d = 6 
condensates. A possible reason for this inconclusive picture may be related to the 
strong correlation between condensates of different dimension. Consider, for example, 
our results for 0\~ A . The 2- and 3-parameter fits lead to very different values since the 
assumptions underlying the two fits are different: in the 2-parameter fit we assumed 
0\q A = 0, whereas the 3-parameter fit preferred the value 0\^ A = 66 GeV 10 and the 
range of values for 0\~ A given in the table is for that fixed central value of 0\^ A . 

It is also interesting to note the agreement of the correlation between Oq~ a and 
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0\- A found in our analysis with corresponding results from [3U [36] • In Ref. [34], the 
linear combination of these two parameters is extracted from weighted finite energy 
sum rules, but no errors were given, while in Ref. Borel sum rules were used to 
find the correlation. In the latter reference, 1-, 2- and 3cr confidence regions for the 
correlations 0\~ A -0\~ A and 0\~ A -0\§ A are presented. The la contours for 0\~ A 
and 0\{~ A are shifted as compared to ours, but the slope agrees well within errors. A 
careful analysis shows that there is also agreement for the 0\~ A -0\^ A correlation with 
the result of Ref. [US]- There, a positive correlation was found from a 2-parameter fit 
which corresponds to fix 0\~ A = 0. With the same assumption we find a correlation 
of the same sign, however a smaller slope. Note that the correlation as shown in Fig. H] 
(left column) appears to be different when fixing Og~ A at its central value which was 
found to be 3.2, i.e. significantly different from zero, in our 3-parameter fit. 

There also exist some previous extractions of QCD condensates in the V and A 
channels, again based on sum rule approaches. The normalisation of spectral functions 
is different from ours and there is also a factor of 8tc 2 absorbed in the definition of the 





Q V-A 


Q V-A 
^8 








-4 ±2.0 


— 1 2 +7 
iZ -n 


78 ±24 


-2.6 ±0.8 


[32]* 


-4.52 ± 1.1 


-10.8 ±6.6 


72 ±28 


-240 ± 100 


m 


J* 


-2.27 ±0.51 


-2.85 ±2.18 


24.1 ±6.1 


-80 ±16 


m 


-8.7 ±2.3 


15.6 ±4.0 


-17.1 ±4.4 


14.7 ±3.7 


m\ 


-7.9 ±1.6 


11.7 ±2.6 


-13.1 ±3.0 


13.2 ±3.3 


m 


-7.2 ± 1.2 


7.8 ±2.5 


-4.4 ±2.8 




[37]* 


-8 ±2 


-2 ±12 






ps 




-6.8 ±2.1 


7±4 






m\ 


-7.7 ±0.8 


11.0 ± 1.0 








-6 ±0.6 


7.5 ± 1.3 






This work 


1-parameter fit [1] 










2-parameter fit [I] 


-6.8lg;° 


o 9+2.8 
°- z -9.2 






3-parameter fit 


o 9+1.6 
°- z -0.4 


-17-0±|;I 


66.0^°;° 





Table 3: Estimated values of the condensates O v d ~ A of dimension d < 12 
in units of 10~ 3 GeV d at leading order. References marked with a * use a 
different normalisation of spectral functions. The values shown are adjusted 
so that they can be compared to those of the present work. 
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A channel 




Of 


^6 


^8 


m 


(1.2... 2.5) 


(4.1... 7.1) 


-(19.1... 13.7) 


[39j 


(1.3... 4.8) 


-(16.5... 1.3) 




This work 


1-parameter fit 


2 fi +0 ' 4 






2-parameter fit 




-1 +a6 

i,u -0.7 




3-parameter fit 


4 8 +L8 


-2 4 +2 -° 


u -°-0.5 



Table 4: Estimated ranges for the dimension d < 8 condensates of the A 
channels in units of 1CT 3 GeV d at leading order. Existing results from the 
literature are presented. Note that the normalisation, for all of them, was 
adjusted so that they can be compared to those from this work. 



condensates. We have translated the results so that they can be compared to ours and 
summarised them for the axial- vector correlator in Tab. HI 

One can remark that the majority of the values found in this work are consistent 
with those from the literature. The sign of Oq, though, disagrees with the vacuum 
saturation approximation and with the results from |llj . 

As a conclusion, we can state that the values and ranges found for the QCD con- 
densates are consistent among themselves and, partly, with previous extractions found 
in the literature even though the agreement between theory and data is very poor in 
the case of the V and V + A channels. Since at present the analyses are still subject 
to a number of restrictions, one can hope that future work will allow us to improve the 
agreement between theory and data further. 

When analysing all four channels, we have assumed chiral symmetry, decoupling 
of heavy quarks and the absence of duality violations. If the chiral symmetry is bro- 
ken, there are also lower-order terms entering the OPE and also mass terms would be 
present both in the OPE and the perturbative expansion. Moreover, there will be also 
a perturbative contribution to the V — A-correlator. Also the inclusion of heavy quarks 
is expected to play an important role at high energies. Their contribution would mod- 
ify the evaluation of the theory prediction Fq CD (s) in Eq. ( TT3l) . It remains to be seen 
whether these effects are negligible or not. 

Duality refers to the assumption that the true function IT(s) can be replaced without 
error by the expression given by the operator product expansion, IIope('S). The term 
duality violation refers to any contribution missed by the substitution IT(s) — > Uope{s). 
As stated already, in our analysis we have assumed that duality violations are absent. 
It is an interesting task to check how the results would change if one would consider 
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duality violating contributions. Unfortunately, little is known about the structure of 
duality violations in QCD and one has to rely on model assumptions like those of Ref. 
[40J. At the time being, possible deviations from duality are suspected to be a major 
source of theoretical uncertainties 141, 421. 
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